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The precise determination of the masses of the neutral and charged heavy Higgs bosons 
is a crucial input for the study of Supersymmetry and its relation with cosmology 
through dark matter. This paper presents a study of e"'"e~ — >■ H'^A'^ and H^H~ pro- 
duction at ^/s=?, TeV. The analysis is performed with full simulation and reconstruction 
accounting for beamstrahlung effects and the overlay of 77 — >■ hadrons events. Results 
I I are presented in terms of the accuracy on the determination of the masses and widths 

^ of the heavy Higgs bosons in two benchmark scenarios. 

H^ 1 Introduction 
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The precise determination of the masses and widths of the neutral and charged heavy Higgs 
bosons is an important part of the study of an extended Higgs sector in models of new 
physics. A non-minimal Higgs sector is one of the simplest extension of the Standard Model 
QQ (SM), realised in the Two Higgs Doublet Model, which acquires a special relevance and 

\^ justification in supersymmetric models. In Supersymmetry the detailed study of the heavy 

^>0 Higgs sector is crucial to assess the relation between new physics and cosmology, through 

^^ dark matter [l]. A high-energy lepton collider is particularly well suited for such a study 

C^ even in those regions of the parameter space where the sensitivity of the LHC becomes 

^^ marginal. In particular, the pair production processes, e^e~ — > H^A^ and e^e~ — > H~^H~ 

^~~^ give access to all four heavy Higgs states almost up to the kinematical limit [2}|4]. 
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2 Simulation 

In this benchmark study we consider two SUSY models with Ma = 742 GeV and 902 GeV, 
respectively, at y^s = 3 TeV and adopt the proposed CLIC beam parameters. The model 
parameters and the masses of the relevant particles are given in Table [Tl The cross sections 
for e+e" — > H'^A'^ and e+e" — )■ H'^H~ pair production are given in Table[2J The properties 
of the heavy Higgs sector are largely independent on the details of supersymmetric models 
and chosen mostly depend on the values of Ma and tan/3. It is interesting to observe that 
the chosen values of these parameters are compatible with the current LHC data and with 
a ~125 GeV lightest supersymmetric Higgs boson, h'^, in the general MSSM. In particular, 
the scenario 1 is compatible with a non-observation of SUSY signals in the channels with 
missing transverse energy by the end of the 8 TeV LHC run and 123< Mh <127 GeV in 
the MSSM [5]. Signal events are generated with Isasugra 7.69 6 and Pythia 6.215 I?]. 
Taking into account the CLIC luminosity spectrum, the predicted signal cross sections for 
H"A° and H+H~^ are respectively 0.7 (0.4) fb and 1.6 (1.1) fb in the scenario 1 (2). The 
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Table 1: Model parameters and Higgs boson masses of the two benchmark points studied. 

Scenario 1 Scenario 2 



Mass Parameters (GeV) 
tan/3, ^0, sign(/z) 
Ma, Mho, Mh± (GeV) 



Ml = 780 M2 = 940 M3 : 
24, -750, + 
902.6, 902.4, 906.3 



540 



mo = 800 mi/2 = 966 

51,0,- 

742.8, 742.0, 747.6 



dominant decay modes are H^ — )■ bh, A^ — )■ bb and i/+ — )■ tb leading to the bbbh and thtb final 
states. The main SM background processes are generated with Pythia 6 . 215 and WHIZARD. 

Table 2: Summary of the processes considered in this study with their production cross 
section and the event generator used. 



Process 


a(fb) 


Generator 


//o^o 


0.7 / 0.4 


ISASUGRA 7.69-HPYTHIA 6.215 


H+H- 


1.6/ 1.1 


ISASUGRA 7.69-HPYTHIA 6.215 


Inclusive 


84.9/ 77.1 


ISASUGRA 7.69+ 


SUSY 




PYTHIA 6.215 


w+w- 


728.2 


PYTHIA 6.215 


z°z° 


54.8 


PYTHIA 6.215 


tt 


30.2 


PYTHIA 6.215 


bbbb 


5.8 


WHIZARD 


w+w-z° 


32.8 


CompHEP+PYTHIA 6.215 


ZOzOzo 


0.5 


CompHEP+PYTHIA 6.215 



In addition, the irreducible inclusive bbbb SM background is generated with CompHep M at 
the parton level and subsequently hadronised with Pythia. We assume the beams to be 
unpolarised. Events are processed through full detector simulation using the Geant-4-based 
Mokka [9] program and reconstructed with Marlin-based 10 processors assuming a version 
of the ILD detector fll], modified for physics at CLIC. Particle tracks are reconstructed by 
the combination of a Time Projection Chamber and a pixellated Si Vertex Tracker. The 
momentum resolution is Sp/p'^ = 2 xlO^^ GeV~^ and the impact parameter resolution 
is (Tij_$ = (2.5 © — (GeV) ) '^^^ '^^'^ parton energy is reconstructed using the Pandora 
particle flow algorithm 12 . Performances for our signal events are discussed in the next 
session. For the 77 — )■ hadrons background simulation two-photon events are generated 
with the Guinea Pig program [13' using the nominal CLIC CDR beam parameters at y/s 
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3 TeV where the hadronic background cross section is modelled following 14 



energies of two colliding photons are stored with their corresponding probability and passed 
to Pythia for the generation of the hadronic events. On average, there are 3.3 77 — ?► hadrons 
events per bunch crossing (BX) with Mjj > 3 GeV. These are passed through the same 
Mokka full detector simulation and overlayed to the particles originating from the primary 
e"'"e~ interaction at the reconstruction stage, using the event overlay feature in the Icio 



persistency package 15 . For this analysis, the 77 background is overlayed only to the signal 



H'^A'^ and H^H events, in order to study its effect on the signal reconstruction and to 
the peaking background from the cross-feed of H^H~ events in the H^A^ analysis and vice 
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versa. 

3 Event Reconstruction and Signal Selection 

The event reconstruction and selection is based on the identification of four heavy parton 
final states in spherical events with large visible energy and equal di-parton invariant masses. 
Details of the bbbb analysis can be found in |4 . Most of the analysis criteria are common to 
both the bbbb and the tbtb channels. The analysis starts with a cut-based event pre-selection. 
Jet clustering is applied to pre-selected events, followed by b and t identification. A kinematic 
fitting is performed to improve the di-jet invariant mass resolution, mitigate the impact of 
machine-induced backgrounds on the parton energy resolution and reject remaining physics 
backgrounds. 

In order to reject particles which are poorly reconstructed or likely to originate from 
underlying 77 events, a set of minimal quality cuts are applied. Only particles with trans- 
verse momentum pt >0.95 GeV are considered, charged particles are also required to have 
at least 12 hits in the tracking detectors and relative momentum uncertainty 6p/p <1. 
The event selection proceeds as follows. First multi-jet hadronic events with large visible 
energy are selected. We require events to have at least 50 charged particles, total recon- 
structed energy exceeding 2.3 TeV, event thrust between 0.62 and 0.91, sphericity larger 
than 0.04 and smaller than 0.75, transverse energy exceeding 1.3 TeV and 3 < N^ds < 5, 
where iVjets is the natural number of jets reconstructed using the Durham clustering algo- 



rithm 16 with 2/cut=0.0025. These cuts remove all the SUSY events with missing energy 



and the e'^e — > // events. For events fulfilling these criteria, we perform the final jet recon- 



struction using the anti-kt clustering algorithm in cylindrical coordinates 17 , implemented 
in the FastJet package, ported as a custom processor in the Marlin analysis framework. 
The choice of cylindrical coordinates is optimal since the 77 — >■ hadrons events are forward 
boosted, similarly to the underlying events in pp collisions at the LHC, for which the anti-kt 
clustering has been conceived and optimised. For each event, the minimum R value at which 
the event has exactly four jets with energies in excess of 150 GeV is used for the clustering. 
The di-jet invariant mass is computed from pairing these jets. Since there are three possible 
permutations for pairing the four energetic jets and the pair-produced bosons are expected 
to be (almost) degenerate in mass, we take the combination minimising the difference AM 
of the two di-jet invariant masses and require |AM| < 160 GeV, 150 GeV for the H'^A'^ 
and H^H^ , respectively. Since the signal events are predominantly produced in the central 
region while the 77 events and most of the SM background processes are forward peaked, 
we only accept events for which the jet with the minimum polar angle, 9, has | cos9\ < 0.92. 
The single most effective event selection to separate signal events from the SM backgrounds 
is ^-tagging, since the signal contains four b hadrons. The irreducible SM bbbb background 
has a cross section of only 0.5 fb and is effectively reduced by the equal di-jet mass constrain 
and kinematic fitting. The 6-tagging procedure is based on the response of the vertexing 
variables of the ZVTOP program flS'. This uses the kinematics and topology of the sec- 
ondary particles in the jet. For this analysis, the vertex variables are supplemented by the 
corresponding kinematic observables for the secondary system reconstructed based on the 
charged particle impact parameters instead of the topological vertexing, when ZVTOP does 
not return any secondary vertex. This procedure increases the efficiency for b jets at the 
higher end of the kinematic spectrum in signal events. Tagging observables are combined 
into a single discriminating variable using the boosted decision tree method with the TMVA 
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package [19^. In the case of charged Higgs bosons, top tagging is also performed. First the 
event is reconstructed as a four jet event and jets are tested for their compatibiHty with 
the top mass. Then a de-clustering procedure is apphed to the jets to study possible jet 
substructure arising from the t — ?> Wb — > qq'b decay. This follows the procedure originally 



developed for identifying highly boosted top quarks at the LHC 20 21 



In order to improve the di-jet mass resolution, we apply a constrained kinematic fit. 
We use the port of the Puf itc kinematic fit algorithm 22 to the Marlin framework. 



Puf itc was originally developed for W^W~ reconstruction in DELPHI at LEP-2 and it has 
been successfully applied for the reconstruction of simulated linear collider events at lower 
energies [3]. The kinematic fit adjusts the momenta of the four jets &spp = apM + bps+cpc, 
where pm is the jet momentum from particle flow, pB and pc are unit vectors transverse 
to Pm and to each other, a, b and c free parameters. For these analyses, we impose the 
constraints Px — Py — 0, E zL \pz\ — \/s and Mjji — Mjj2, where the second condition 
accounts for beamstrahlung photons radiated along the beam axis. We only accept events 
with a kinematic fit x^ < 5. After the kinematic fit, the relative jet energy resolution 
RMSgo/^jct for 6 jets is 0.075±0.002 and 0.096±0.002 without and with 77 background 
overlayed, respectively. The use of a kinematic fit also mitigates the effect of the overlayed 
77 events on the di-jet mass resolution. Since we do impose the nominal centre-of-mass 
energy, allowing for beamstrahlung, jet energies are rescaled in the fit to be consistent with 
that of -^/s. The di-jet invariant mass resolution changes from (38±4) GeV, for the raw 
particle flow objects, to (29±3) GeV, using the kinematic fit, and to (18±2) GeV, imposing 
the equal mass constrain, in absence of 77 background. Using the kinematic fit and equal 
mass constrain we obtain (19.3±3.0) GeV with the 77 background overlayed. The use of the 
equal mass constrain is justified in the case of the e+e" — > H'^A'^ channel in SUSY models. 
In fact, the mass splitting between the A'^ and H'^ bosons is always smaller than both the 
natural widths and the experimental di-jet mass resolution, in particular in the region of 
MSSM parameters where their determination is most crucial to fix the relic dark matter 

The 
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density Q.^h . This has been verified through a dedicated study using FeynHiggs 

few cases in which the mass splitting becomes large are also characterised by very large 

natural widths of the two bosons. 



4 Results 

The di-jet invariant mass distributions for the bbbb and tbtb final states are shown in Figurelll 
The masses and widths of the heavy bosons are extracted from a multi-parameter fit to 

Table 3: Summary of the fit results for the masses and widths for benchmark scenario 1, 
without and with 77 background overlay. The quoted uncertainties are statistical only. 





State 


Mass 


Width 






(GeV) 


(GeV) 


No 77 


A7iJ" 


742.7±1.4 


21.7±3.3 


No 77 


iJ± 


744.3±2.0 


17.0±4.7 


77 (20 BX) 


A^H" 


743.7±1.7 


22.2±3.8 


77 (20 BX) 


H± 


746.9±2.1 


21.4±4.9 
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Table 4: Summary of the fit results for the masses and widths for benchmark scenario 2, 
without and with 77 background overlay. The quoted uncertainties are statistical only. 





State 


Mass 


Width 






(GeV) 


(GeV) 


No 77 


A^/H" 


902.1±i.9 


21.4±5.0 


No 77 


iJ± 


901.4±i.9 


i8.9±4.4 


77 (20 BX) 


A7i/o 


904.5±2.8 


20.6±6.3 


77 (20 BX) 


H^ 


902.6±2.4 


20.2±5.4 
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Figure 1: Di-jet invariant mass distribution of fully simulated and reconstructed events in 

the channel e'^e^ —>■ bbbb (upper row) and e^e^ — > tbtb (lower row). The contributions from 

the different processes are shown. The dotted line shows the fitted background distribution 

which includes the peaking contribution of the cross feeds from the i/'^j4° and H^H~ 

processes. The continuous line shows the global fit including the signal. 



these distributions. In order to assess the effect of the 77 — > hadrons background on the 
accuracy of the mass and width reconstruction, we repeat the analysis assuming first no 77 
background and then we overlay the background assuming an effective 10 ns time stamping 
capability of the detector, corresponding to the integration of 20 bunch crossings (BX). 
Results are summarised in Tables |3] and lU 

We observe that the heavy Higgs mass can be measured to a relative statistical accuracy 
of better than 0.5%. This accuracy is preserved in the presence of 77 background by apply- 
ing the anti-kt jet clustering and kinematic fitting, provided the global detector stamping 
resolution is <10 ns. The width can be measured to a relative accuracy of 0.20-0.30. These 
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measurements are important both for constraining the model parameters through Ma and 
tan/3 and for assessing the contribution of the x°x° -^ A'^ ^ bb pole annihilation process 
in setting the relic dark matter density in the Universe [I]. These accuracies ensures that 
the heavy Higgs mass contribution to the statistical precision in the extraction of fl^h^ in a 
generic MSSM model, with dark matter annihilation through the A'^ pole, is of order of 0.10. 
In the specific case of scenario 2, where tan/3 is large and Af^=742 GeV, these constraints 
on r^o and r^± also yield a relative statistical accuracy on the determination of tan/3 to 
-0.20. 

5 Conclusions 

A study of H''^A'^ and H^H^ production in y^ — 3 TeV e+e^ collisions performed on 
fully simulated and reconstructed events with realistic backgrounds shows that the mass 
and widths of these bosons can be determined with a relative accuracy of 0.05 and 0.2- 
0.30, respectively, provided time stamping with a resolution <10 ns can be achieved in the 
detector. 
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